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In order to describe the thermodynamics of the glassy systems it has been recently introduced 
an extra parameter, the effective temperature which generalizes the fluctuation-dissipation theorem 
^^ ' (FDT) to off-equilibrium systems and supposedly describes thermal fluctuations around the aging 

state. Using this concept we investigate the applicability of a zeroth thermodynamic law for non- 
equilibrium systems. In particular we study two coupled systems of harmonic oscillators with Monte 
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JH ' Carlo dynamics. We analyze in detail two types of dynamics: 1) sequential dynamics where the 

coupling between the subsystems comes only from the Hamiltonian and 2) parallel dynamics where 
there is a further coupling between the subsystems arising from the dynamics. We show that the 
coupling described in the first case is not enough to make asymptotically the effective temperatures 
of the two interacting subsystems equalize, the reason being the too small thermal conductivity 
between them in the aging state. This explains why different interacting degrees of freedom in 

M ' structural glasses may stay at different effective temperatures without never mutually thermalizing. 

Ch ' 

I. INTRODUCTION 

C^ i The dynamics of glassy systems has been a subject of intensive researchB. Despite the fact that glassy systems are 

ofF-equilibrium systems, some regularities that allow the rationalization of the problem have been found. One of the 
most striking regularities is the presence of aging. This means that the correlation and rjesponse functions are not 

S^ only functions of time-differences but also of the time elapsed since the system was prepared^. Thus, qualitatively, the 
longer one waits in the low temperature phase, the smaller the response to an external field will be. A salient feature 
of systems in equilibrium is the fact that the linear response functions and the equilibrium fluctuations are related 
by the well known fluctuation-dissipation theorem (FDT)EJ. This relation does not hold for off-equilibrium systems. 
Several studies of spin-glass mean-field models have shown that a generalizatioiiiDf the fiuctuation-dissipation theorem 
is possible through the definition of the "fluctuation-dissipation ratio" (FDR)Ij'Q: 
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§: ^(M) = ^|^(t>.), (1) 

^^ which is equal to 1 in equilibrium. It turns out that the behavior of the quantity X{t^s) is non trivial in the 

2^ limit t,s ^ oo. If the lowest time s is sent to inflnity the quantity X(t,s) becomes a non-trivial function of the 

■"".^ ] autocojxelation C{t, s). This a strong statement which has been proved to hold in the framework of mean- field spin 

jrt ■ glassegju. Moreover, it has been recently recognized that the quantity X is generally related to the Parisi order 

Cj ' parameter P{q) whick appears in equilibrium studies of spin-glasses providing a natural link between the static and 

I , dynamical propertiesQ. 

'^ ' What is the physical interpretation of X7 According to relation (|l|) the fluctuation-dissipation relation would be 

Ch , satisfied if the temperature into the_cight hand side of (||) were T/X{t, s). This last ratio receives the name of effective 

O ■ temperature and it has been shownO that it has some of the good properties of a macroscopic temperature. In fact a 

proper thermometer coupled to the slow degrees of freedom can measure it. The value of Tcs(t, s) = T/X{t, s) would 

then be different (and higher) than that of the thermal bath. The question about the convenience of this temperature 

k>( ■ to describe the non-equilibrium behavior has been ajSiibject of controversy in the last yearall. While there are some 

''T^ . evidences (not only theoretical but also experimentaOu) that the violation of FDT gives a good temperature in the 

d ' thermodynamic sense, it is unclear what properties of standard (i.e. equilibrium) temperatures are common to the 

non-equilibrium ones. 

The motivation of this paper is to answer to the following question: How effective temperatures equalize when 
two systems out of equilibrium are put in contact? In other words, does there exist a zeroth law for non-equilibrium 
systems? Let us imagine about a vitrified piece of silica quenched to the room temperature. Because the glass is 
off-equilibrium its effective temperature is higher than room temperature. But, if we touch the piece of glass it is 
not hotter than the room temperature. We must conclude that some degrees of freedom within the piece of silica 
are thermalized to the room temperature while other remain non-thermalized and still hotter. Touching the piece of 
silica we feel the fast modes, not the slow ones. This poses the question, how is that possible that different interacting 
degrees of freedom have not reached thermal equilibrium for sufficient long times? Despite of some considerations 



present in the literatureOll3 there are no clear answers to this question. We beheve that some of them may require 
a more deep understanding through a detailed analysis of an illustrative example as a previous stage to offer more 
simple and generic considerations. It is our purpose here to follow this route trying to give a general answer to this 
question by deriving exact results in the framework of a solvable model. l-l 

The model is a set of harmonic oscillators evolving by Monte Carlo dynamics introduced inS (hereafter referred as 
BPR model). The importance of this model relies on the fact that it is exactly solvable and shows one of the main 
features of glasses, namely aging in correlation and response functions. Our interest will be in considering two coupled 
sets of harmonic oscillators. Thus, we can see how the main observables are affected by the coupling, in particular how 
the effective temperature evolves for the two sets of interacting degrees of freedom (represented by the two different 
sets of harmonic oscillators). The interaction may then appear through the Hamiltonian or through the Monte Carlo 
dynamics itself. We will discover that the effective temperature for the two sets of oscillators depends on how the 
coupling is done, and we will understand why in vitreous systems different degrees of freedom may stay at different 
temperatures without thermalising at very long times. The central idea is that interacting non-equilibrium systems 
each one with very different effective temperatures may not equalize because the conductivity in the aging state can 
be extremely small. In this sense the utility of the extension of the zeroth thermodynamic law to the non-equilibrium 
aging state is questioned due to the smallness of the non-equilibrium conductivities. 

The paper is organized as follows. Section II describes the main aspects as well as the interest of the model. Section 
III describes the two classes of couplings we have considered. Section IV analyzes the case in which the main coupling 
is ruled by the Monte Carlo dynamics. Section V describes the case where coupling appears only in the Hamiltonian. 
Sections IV and V show how to solve the dynamics of the system. The reader who is not interested in technical issues 
can skip them. Section VI discusses the results and the physical consequences of our work. The last section presents 
the conclusions. Three appendices are devoted to some other technical issues. 

II. A SIMPLE AND SOLVABLE MODEL OF GLASS 

As a simple model of glass we will consider a system of uncoupled harmonic oscillators evolving with Monte Carlo 
dynamics. The Hamiltonian is: 



\kj:x^ . (2) 
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This model was introduced incl and was also reviewed iimU. 

The low-temperature Monte Carlo dynamics of an ensemble of linear harmonic oscillators shows typical non- 
equilibrium features of glassy systems like aging in the correlation and response functions. The interest of this 
model is that the slow dynamics at low temperatures is a consequence of the entropy barriers generated by the low ac- 
ceptance rate. The simplicity of this model makes it exactly solvable yielding a lot of results about the non-equilibrium 
behavior. 

The Monte Carlo move consists on the following: the Xi are moved to Xi + ri/\/N where r^ are random variables 
Gaussian distributed with zero average and variance A^. The move is accepted according to the transition probability 
W{/S.E) which satisfies detailed balance: W{ls.E) = W^(— A£') exp(— /3Ai<^), where A£' is the change in the Hamil- 
tonian. In Appendix A we show the computation of the correlation and response functions. Here we only quote the 
main results, 

1. Slow decay of the energy. The evolution equation for the energy is Markovian. This simplicity allows for 
an asymptotic large-time expansion showing that the energy decays logarithmically E{t) ^ l/log(i) and the 
acceptance ratio decays faster A{t) ~ l/{t log(i)). 

2. Aging in correlations and responses. The correlation function C(i, s) is defined by: 
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^(*'*) = ]^E^'W^*w • (3) 

The response function is calculated by applying an external field to the system. Then, the response function is 
the variation of the magnetization of the system when the field is applied: 



with the magnetization given by, 

N 



Details o n how to solve correlations and responses are given in Appendix A. The final results are equations 



(102,111). Both correlation and responses show dominant s/t scaling with logarithmic corrections. The asymp- 



totic scaling behavior is given by. 



C{t,s)^C{s,s)^ , G{t,s)=G{s,s)^e{t-s) , (6) 

with C{s,s) = ^^' , G{s, s) ~ ^-^ where the expression f{t) is given in equation (98) and L(t) ^ i(log'^(t))ta. 
The slow decay of the response £inction shows the presence of long-term memory which manifests as aging in 
the integrated response functions. 

3. The effective temperature. As said in the introduction, the effective temperature is defined in terms of the FDR 
eq-i): 

dC(t,s) 

G{t,s) 

In equilibrium E{s) — T/2 and we recover the expected result Toff — T. Interestingly (R) yields a result for 
Teff which only depends on the smallest time s. The unique dependence of the effective temperature on the 
lowest time s is generally believed to be satisfied in the asymptotic large s limit for generic structural glasses 
and spin-glass models with a one step of replica symmetry breaking. This expectation holds here for all times. 
At zero temperature when slow motion sets in, the system never reaches the ground state and ages forever. In 
this regime the effective temperature verifies in the long-time limit (i.e s — > oo ): 

T^s) = 2Eis) + y^^ ^ ms), (8) 

This gives a thermodynamic relationship between the effective temperature and the dynamical energy in the off- 
equilibrium regime showing how the equipartition theorem can be extended to the glassy regime. The effective 
temperature measures how a quasi-stationary or adiabatic hypothesis is exact for the present model suggesting 
that some features of equilibrium thermodynamics may be applied to the aging regime. 

III. TWO COUPLED SYSTEMS 

Now we consider the case in which we couple two systems of harmonic oscillators. In this case it is possible to 
compute analytically how one system affects the other without loosing the benefit of evaluating Gaussian integrals. 
The Hamiltonian we have to deal with is: 

^ = fE-? + fE^.^-^E-»^^' (9) 

4=1 i=l 1=1 

where we take K1K2 > e^, otherwise the system has no bounded ground state. We define the following extensive 
quantities (per oscillator): 

i—l i—l i—1 

where Ei and E2 are the energy of the bare systems while Q is the overlap between them. In this case we also 
consider Monte-Carlo dynamics, where the transition probability is performed by the Metropolis algorithm which 
satisfies detailed balance. The random changes in the degrees of freedom cc^, yi are defined in the same way we have 
explained in the previous section for the case of a single system. But there are different ways to implement the 
dynamics in the model depending on the updating procedure of the variables Xi,yi. Here we have analyzed two 
important and different procedures which yield quite different results: 



1. Uncoupled or sequential dynamics. In this case the two sets of variables x and y are sequentially updated. First 
the Xi variables are updated and the move is accepted according to the total change of energy AE — AEi — eAQ. 
Next, the variables yi are changed and the move accepted according to the energy change AE — AE2 — eAQ. 
This procedure is then iterated. In this case, the dynamics does not affect simultaneously the two sets of variables 
but each set is updated independently from the other. The only coupling between the two sets of oscillators 
comes from the explicit coupling term eQ in the Hamiltonian. Note that for e = the dynamics becomes trivial 
because the dynamical evolutions are that of two independent sets of harmonic oscillators everything reducing 
to the original model described in section II. 

2. Coupled or parallel dynamics. In this first case the Xi,yi variables are updated in parallel according to the 
rule Xi ^ Xi + ri/\/N, yi ^ yi + Si/\/N. The transition probability for that move W{AE) is determined by 
the change in the total energy AE — AEi + AE2 — eAQ introducing, on top of the explicit coupling term eQ 
in the Hamiltonian, an additional coupling between the whole set of oscillators through the parallel updating 
dynamics. Contrarily to the uncoupled case, the e = case is interesting by itself because it shows how this 
kind of dynamical coupling strongly influences the glassy behavior. In fact, in the limiting case e = 0, there 
will be some changes which make the energy of one of the two systems increase, this change being accepted 
because the total energy will decrease. Because of that, despite of the fact that there is no direct coupling in 
the Hamiltonian the dynamics turns out to be strongly coupled. 

In what follows we describe the main set of quantities we are interested in. The solution of the dynamical equations 
for the coupled and uncoupled cases is very similar. The Appendix B shows in detail the derivation of the dynamical 
solution for the uncoupled case. 

A. Correlation, overlaps and responses 

On top of the time evolution of one-time quantities our interest will also focus on the behavior of two-times 
quantities such as correlations and responses. These quantities will refer to three classes of systems: the set of 
oscillators described by the x variables, the set of oscillators described by the y variables and the whole set of x and y 
variables. In the rest of the paper, as a rule, the subindex 1 will refer to quantities describing the set x of oscillators, 
the subindex 2 will refer to quantities describing the set y of oscillators and the subindex T will refer to quantities 
describing the whole set of oscillators x plus y. The main set of correlation and response functions we are interested 
in are: 

• Correlations. The correlation function for the sets x and y, 

1 ^ 1 ^ 

^i(^'^) = ]^E^^W^^W ' C'2(t,s) = -^y.(t)2/,(s) , (11) 

as well as the global correlation Crit, s) — ^(Ci(i, s) + C2{t, s)). 

• Overlaps. These are cross-correlations involving different sets of variables: 

1 ^ 1 ^ 

^i(*'^) = ]^IIy*W^*(^) ' Q2(i,s) = ^^x,(i)y,(s) , (12) 

with Qi{t, s) = Q2{s, t). As we will see later, it is useful to define these two functions Qi, Q2 which essentially 
are the same overlap function but acting on different time sectors. 

• Response functions. The response function for the sets x and y are defined in the following way. Define the 
magnetizations for the two sets of oscillators x and y, 

1 ^ 1 ^ 

i=l i=l 

Consider also two external fields hi and ft,2 conjugated respectively to Mi and M2, 



H ^ ^^xi + ^^y'i-^{hiXi + h2yi)-e^x,yi . (14) 

z— 1 i— 1 i i 

We define four types of response functions 6*1^2, G[ 2- The Gi{t, s), G2{t, s) functions measure the change in the 
magnetization Mi{t), M2{t) induced by their respective conjugated field /ii,/i2 applied at time s. These are 
defined by 

-<-'^(^),„. • 

where the index i = 1, 2 represents each one of the systems. Apart from these two response functions we may 
define the global response function GT{t,s) as the change in the global magnetization Mt = ^(-^^1 + M2) 
induced by a field conjugate to the total magnetization, 

The primed response functions G'i{t, s), G2(t, s) functions measure the change in the magnetization in each set of 
oscillators Mi(t), M2{t) induced by a conjugated field (respectively h2,hi) apphed on the other set of oscillators 
at a time s: 

where the indices i,j = 1,2 are different i 7^ j. In the absence of a coupling term eQ in the Hamiltonian (|9|) the 
two response functions G[ 2 vanish but for e 7^ they enter into the solution of the dynamical equations. 

• Effective temperatures. From the correlation and response functions we may define three effective temperatures: 
T^ff for the system 1, T^g for system 2 and T^ff for the global system. These are defined as follows, 

/ dCl{t,s} \ / dC2{t,s) \ / dCT{t,s) \ 

We will analyze in detail the three effective temperatures for the coupled and the uncoupled cases. From them 
we will learn whether the systems equalize their temperatures and how they do. 

B. Equilibrium regime 

Here we present the results for the statics for the general model (0). The equilibrium solution is the stationary 
state of the dynamics coinciding for both coupled and uncoupled dynamics. The results for the one-time quantities 
can be simply evaluated from the partition function, 

/oo 
dxdyexp{-/3H), (19) 
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which involves simple Gaussian integration. By performing the appropriate partial derivatives we calculate the 
different thermodynamic quantities: 

£,eq _ ^eq _ K,K2T ne^-'l- Ml! (20) 

where the parameters /, J are defined by, 

/ = — !-A^ + ^-A? , J=— !-^ + ^-^ , (21) 
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where the total energy E = Ei + E2 — eQ is given by the equipartition relation E ^ T. Note that K1K2 — e^ > in 
order for E'^, E'^ to be positive. 

The equilibrium correlations C, overlaps Q and responses G, G' only depend on the time differences. While the 
precise form of these functions depends on the particular type of dynamics, the magnetic susceptibilities do not. These 
are given by: 

/■OO T^ pOO J^ POO -t 

XI = / Gi(t)di - ^ ^ , X2 = / G^mt - ^ , , XT = \ GT{t)dt = -(xi + X2), (22) 

Jo KxK2 - e-^ Jo KxK2 - e^ Jo 2 

and are temperature independent as expected for oscillator systems. Nonetheless, in equilibrium the three effective 
temperatures (^_8|) coincide with the bath temperature T . 

IV. THE DYNAMICALLY UNCOUPLED (OR SEQUENTIAL) CASE 

In this section we solve the dynamics of the thermodynamic relevant quantities for the case in which the two 
subsystems of oscillators are dynamically uncoupled. As explained in the previous section, in this case we make a 
sequential dynamics avoiding direct dynamical coupling effects coming from the Monte Carlo dynamics. The derivation 
of the dy namical equa tions is explained in the Appendix B. The equations for the energies and overlap ([lO| ) are written 
down in ( p^ , p^p3| ), 

^^' <^E, - eQ)/«,(t) + 1 (l^ + "^ercfiaA (23) 



dt ' " -^/^mv > • 2\ [3 

^ = H2E2 - eQ)fnAt) + \ (^ + ^erc/(a.)) (24) 

with the fohowing definitions: 

Ri^Ei^eQ+—^ , R2^E2^eQ + -—^, (26) 

A1A2 J^1^2 



/fl.(i) = ^^/3exp('-/3^^(l-2i?,(i)/3)')er/c(a,(i)(4i?,(i)/3-l)) with a, ^ J ^^ , (27) 

where the error function was defined in (|97|). 

Definitions ( p7|) hold for i = 1,2, each i representing one of the two systems. Note that the whole dynamics is 
contained in the function /ij.(t) . In what follows we will be especially interested in the zero-temperature case where 
relaxation time diverges and dynamics is slow and glassy. For T = the function fR-(t) in ( |27| ) becomes 

/^. = ^exp(-a2) . (28) 



A. Asymptotic long-time expansion for the one-time quantities 



The asymptotic solution of equations (^3|J24|j25|) may be guessed from the behavior of the energy (113) for the 
uncoupled systems. Trying a solution of the type 

log(t) log(i) log(t) 

we can only solve the asymptotic behavior in the limit e « 0. This is a consequence of the fact that the quantities 
i?i and i?2 are different in general and we have a system of four equations with three parameters. There is not any 
general solution for this system, but in the limit e ~ the quantities i?i and R2 become identical to the corresponding 
energies to leading order yielding only three equations with three unknown parameters (a, 6, c). In this limit the value 
of the coefficients a, 6, c may be easily obtained yielding 
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where J was defined in eq.(pi|).At first order in logaritlimic corrections 1/ log(i) we find in the limit e w : 
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16{KiK2 - e2) log(i) 



16iKiK2 - e2) log(i) 
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eJ 



8{KiK2 - £2) log(i) 



(30) 



(31) 



Note that, in the long time limit both energies Ei,E2 tend to zero logarithmically but their relative difference 
^^ ^ stays finite. In the limit of small coupling constant we can do a more refined expansion yielding: 



E2 = 



KIK2AI 
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l&{KiK2 - £2) log(t) + \ l0g(l0g(i)) 



KiKjAl 1 



16(i^ii^2-e2)log(i) + ilog(log(i)) 
eJ 1 

^8(ifii^2-e2)log(t) + llog(log(t)) 
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(32) 
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where we have put explicitly the terms of order 1/ log^(i) as sub-dominant corrections. These terms come from the fact 
that the true expressions for the energies and the overlap should be, in order to match the coefficients in (E3,E4,G5[): 



Ei = 



\og{At) 



E2 = 



logiBt) 



Q = 



log{Ct)' 



(35) 



which gives a correction of order l/log^(i) in expressions (|3l|). 

In Fig.n^ we show the evolution for the energies and the overlap for two systems of harmonic oscillators with a small 
value of e. We also show the asymptotic behavior (|3^j3^,^). We can see that the two energies remain different even 
at long times. We will see that this feature is very important for describing the non-equilibrium state of the whole 
system. We can also see that the asymptotic expansions are in good agreement with the numerical solution of the 
dynamic equations. Nevertheless, there are systematic deviations at long times being consequence of the limited range 
of validity (e <C 1) of the asymptotic solution ( |3^ , |33| , |34| ) . If KiA^ = K2A2 the energies of the two oscillators become 
identical (note that equations (£^,^,Q) only depend on the constant KA^) and asymptotic dynamics also. 
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FIG. 1. Tlie decay of tfie energies and tfie overlap for two systems with K\ = 2, K2 = 1, Ai = 1, A2 = 1 and e = 0.2. Tlie 
longest lines are the numerical solution for the dynamic equations, while the shorter ones are the corresponding asymptotic 
behaviors. 



B. Correlations and responses. 

The set of equations for the four correlation and overlap functions defined in (^,^2|) can be written as: 
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(36) 



with the subsidiary boundary conditions 

2^1 (s) 



Ciis,s) 
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C2(S,S) 



2^2 (g) 



Qi{s, s) = Q2(s, s) = (3(s) 



(37) 



The equihbrium solution can be appropriately worked out because the matrix coefficients are time-independent. If 
we write the matrix equation in compact form 



dt 



MC 



with C = (Ci, C2, Qi, Q2) the solution is 

C(t,s) = C(s,s)exp(Meq(i-s)) 



(38) 



(39) 



The precise results for correlations and overlaps are reported in the Appendix B (formulae ( |151 -159)), the initial 
conditions being given in (j^Tjj. In the non-equilibrium case it is not possible to write down an exact solution for 
equation (pSl) for any value of e. The formal solution of (pSl) is 



(40) 



C(i,s) =rexp(e / Bi{s'))C{s 

J s 



which can be worked out perturbatively up to any order around e = (T stands for the time ordered product). In the 
Appendix C, we give some details how to construct such expansion. Up to order e^ the solution for the components 
Ci , C2 of the four component vector C are 

Ci(i,s)=exp(^- / fR,{x)dxj[^ ^^ ^'^~K^ dt fR,{t)eJ^ '''''' ' '"'' " + 

df'fR, {t')fn. (i")e^*'' ^^"^ ^^^^^"^ ^^^^'^) , 



K2KI 



(41) 



C2[t. s) = exp(- f^ fRA-)dx) (^ -f e9^ f^ dt'fn,it')ei:^f^^^^^~f'^^^^-^^'^ + 

e' ^^ /* dt' f df'fR, {t')fn, {t")eC W«^ f^)-^-. (-))'^^) . (42) 

Similar expansions are obtained from Qi{t, s), Q2{t, s). Here we do not report them because correlations are enough 
to analyze the effective temperatures. Similarly we can also obtained expressions for the responses as detailed in the 
Appendix B. The time evolution for the four possible response functions is given by 



^^^--{G,{t,s)fnAt)-^Sit-s)-^^G',{t,s)fnAt) 
^^§^ = - (G2{t,s}fR,{t) - l^S{t - .) - S.G[{t,s)fnAt)) 
^^^ - - (G\it,s)fnAt) - -^G2{t.s)fnM 
^^^ = - (G',{t,s)hAt) S.G,it,s)fnJt) 



(43) 
(44) 
(45) 
(46) 



As explained in Appendix B these equations must be solved with the subsidiary boundary conditions, 



G^is,s) = i^ ; G2(s,s) = ^^ , Gl(s,s) = , G',{s,s)=0 . (47) 

Ai A2 

The initial conditions for G'i,C?2 come from the deha-terms in their equations. The other two initial conditions for 
G'l, G'2 come from the fact that there is no discontinuous jump in the response function of one system when we apply 
the field to the other system. This result also holds in the framework of the Langevin dynamics and manifests in the 
equations for the magnetizations (see in Appendix B (143,143)) as the absence of a field /12 in the equation for Mi 



and the absence of a term hi in the equation for M2. 

In equilibrium the expressions for the responses Gi, G2 are given in the Appendix B. Up to order e^ the expression 
for the off-equilibrium responses Gi , G2 can be solved analogously as done for the correlations and are given by 

Gi{t,s)^^eMffRA^)dx){l + -f-fdt' f dt'7«,(t')//..(i")e^'''^^"^'^^-^"^^^^^'") (48) 

G2 {t, s) = ^ exp(^' fn, {x)dx) (l + ^ j' dt' j' dt" Jn. {t')fn, (t")e/*''' ^^"^ ^^^"^"^ (^»'^) . (49) 



V. THE DYNAMICALLY COUPLED (OR PARALLEL) CASE 

For the dynamically coupled case the calculations proceed similarly as to the previous dynamically uncoupled case. 
The evolution equations for the overlap and the energies (M) are: 






dt \Q-j)fEAt) (50) 



dt 4 -^^Tmj^y-^-'''^^^'''^'''^^^^)-^ ^''^ 

dE2 K^Aj^^.J J \ fK^Al ,^ ^2.^ , ,^A2.. \ /s«(t) 



3t 4 ^'^^^iTii^jni^-'^^^^^^^^^^^^^j^^T^ • ^''^ 

I and J were defined in (p^) and the new quantities E^q, /^^ and Et {Et is not the total energy) are given by 

En =Et~ 2eQJ + e^I , Et = EiKiAJ + E2K2AI (53) 

/b„ {t) = J/3 exp {P^En - J/3) erfc ( ^^"^'^ ^ (54) 
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and the expression (54) at zero-temperature is 



/^„(t) = ^^exp(--^) . (55) 



A. Asymptotic long-time expansion for the one-time quantities 

Proceeding similarly as done in the former section we can find the asymptotic expressions for the energies and 
overlaps. In this case we can find a solution for e finite due to the fact that in this case we have only one dynamic 
function /^^ . We find, to leading order in 1/ log(t) 

^ K1K2J ^ K1K2J ^ eJ 

' 8{KiK2-e^)log{t) ' ' 8{KiK2-e^)log{t) ' ^ ^1X2 ~ e^)log{t) ' ^°"^ 

Note that, contrarily to results (pi|) for the dynamically uncoupled case, the energies Ei, E2 asymptotically coincide 
and the relative difference {Ei — E2)/Ei vanishes in the long-time limit. This difference of behaviors is not casual 
and has a physical interpretation that we will discuss later. The more precise expansion turns out to be. 



E, = 



Eo^ 



Q 



K1K2J 1 

8(ifiif2-e2)log(t) + ilog(log(i)) 



K1K2J 



1 



0{ 
0{ 



1 

log'(i) 
1 



8(iflif2 - £2) log(t) + i l0g(l0g(t)) ■ - '\0g\t) 



) 



4(Xiir2-e2)log(t) + ilog(log(t)) 



^log^W 



(57) 
(58) 
(59) 
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FIG. 2. The decay of the energies and the overlap for two systems with Ki — 2, K2 ~ 1, Ai — 1, A2 = 1 and e — 0.3. Black 
lines are the numerical solution for the dynamic equations, while the blue ones are the corresponding asymptotic behaviors. 



Let us stress that, contrarily to the dynamicaUy uncoupled or sequential case the previous expressions are valid 
to any order in e. The origin of the 1/log (t) terms in previous expressions is the same as in the uncoupled case 
(^,^|4|). In Fig. H, we show the numerical solution for the evolution of the energies and the overlap as well as the 
asymptotic expansions ( |57| , |5q , p^ . 

We have said that the relative difference {Ei — E2)/Ei vanishes in the long-time limit. It is not difficult to see 
how this happens. The time-evolution for the quantity E1/E2 is easy to derive from eqs.(^,^,|2|) in the asymptotic 
long-time limit En — > 0. One then finds the following expansion to leading order 



^.1 
E2 



2Jlog(i) 



(60) 



K2A2 the correction is even smaller. 



If KiAj 

Interestingly, this leading correction does not depend on e showing that the two energies Ei , E2 approach each 
other at a rate determined by the fact that the whole dynamics of the model is coupled and not by the fact that the 
two oscillator systems are coupled by the presence of a term eQ in the Hamiltonian. The behavior of this quantity is 
shown in figure (H) for different values of e together with the asymptotic expansion (RQ) . 



e = 0.1 -0.9 
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time 



10 



10 



10 



FIG. 3. Relative energy diflterence Ex/E2-l for two systems with K^ = 2, 7^2 = 1, Ai = 1, A2 = 1 and e = 0.1, 0.3, 0.5, 0.7, 0.9 
(from top to bottom). The asymptotic prediction (pOJ) is also shown. 



B. Correlations and responses. 

Following similar methods as for the dynamically uncoupled case presented before we can write down the equations 
for correlations and overlaps 

^^^ = - {K,^lC,{t,s) - eA?Qi(M)) %^ (61) 

^^ = - (A^.A^C.(M) - eA^Q.(M)) ^ (62) 

^^^ = - (i^2A^Qi(t,.) - .A^Ci(t,.)) l^ (63) 

^^ = - (i^i A?g.(t, .) - .A?C.(i, .)) %^ , (64) 

with the subsidiary boundary conditions given by equations (|37| ). 

In matrix form these equations reduce to the equation (|38|). As explained for the uncoupled case, this set of 
equations can be exactly solved only in the equilibrium regime where the coefficients are time-independent. The 



solution is then given by the equation (|39|), the expressions for correlations and overlaps are the formulae (151 156) 



KlAf , , _ eAf . . „ _ eA^ . , _ K^C^i 



K2Ai 



with a = ^fEM\h = -jffEMiC = -2ffEM',d = tj^fEM and the formulae ( [152| , |157D with a ^ fj^/fiiv; ^ = 

eA| p . _ tAf r 7 _ KiAJ r 

2J JEn ' ^ 2 J -fEN ' '^ 2 J J En • 

In the most general case where the coefficients of the matrix equation are time dependent the exact solution can 
be written in the closed form (^) which can be expanded to any order in e as explained in the Appendix C. As in 
the uncoupled case we present here the results up to order e^ only for the correlations. 



Responses Gi^2j G'l 2 can be worked out in a similar way as shown in the Appendix B for the uncoupled case, 
dGi{t,s) _ KiAl 



Ci (i, s) = exp 

V ^-J Js 

C2{t,s) = exp 



(^Gi (i, s)fEN (t) - ^^S{t -s)- -^G'2 {t, s)fEN (i)) (67) 

(68) 



dt 2J 

'-^ = -^{G2it,s)fEjt)^^Sit-s)-^^G[it.s)fEAt)) 

^^T^ = -^ (G[{t,s)fEAt) - -|-G2(i,.)/^„(i)) (69) 

^^m^ = -^ (G'2(M)/^„(t) - ^^G,it,.s)fEAt)) , (70) 

with the subsidiary boundary conditions (note that for Gi , G2 they are different from those in (|4 
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Gi{s,s)^^I^ , G2{s,s) = ^^ , G[{s,s)^0 , G',{s,s)^0 . (71) 

It is a simple exercise to check in equilibrium whether these responses give the correct value of the susceptibility 
(P2|). In equilibrium responses only depend on the difference of times. As the susceptibility is just the integral of the 
response function we can integrate the equations (for simplicity we shall consider that s = 0). Then, the equilibrium 
susceptibility of every system is just: 

/•OO 

X = / G{t)dt . (72) 

Jo 

Integrating the equations for response functions we obtain: 

Gi(oo) - cm = ^j— [xi + ^X.j = ~^J]^ (73) 

Gi(oo) - G'M = -^^^^ (x'. + ^Xi) = (74) 

c.(^)-c.(o)^-^^^(..^x;)^-^^ (-) 

G[i^) - G'M = - ^%^^^' (x[ + ^X^) = . (76) 

At very long times, ergodicity imposes Gi,2(oo) — G'i2{°°) — 0- These equations give the exact results for the 
equilibrium susceptibilities eqs.(p2). In equilibrium the responses can be easily computed and one gets (to keep 
formulae at minimum we only report the results for Gi and G2), 

Gi(t, .) = ^%^ (^-^ exp(A,(t - .)) - A_iL exp(A2(t - .))) (77) 

G2(M) = ^%^(^exp(Ai(i-.))-^exp(A2(t-.))) , (78) 



with the usual expression ( |155[ ) for Ai, A2. For Gi the values of the constants are: 

« = 2J~ "' "2j" "' '^^"2j" "' 2J~ " C^^-* 



while for G2 the same results (|79|) are valid but interchanging the indices 1 and 2. 

In the general off-equilibrium case the result for Gi , G2 , G'l , G2 can be worked out perturbatively. Here we only 
write the expression up to order e^ 



2 / &-. A 2 /•* 



1+6^ 



Gi(t,s) = ^^exp(-^y^ /^„(x)dx) X 
^j^dt'j^ dt"fEAt')fE.{t")e^ ^ )/../-«(-)'^-) (80) 

G2(M) = ^-P(-^ f fEA^)dx) X 



(l + e^^y dt'y di"/i.„(t')/i.«(i")e( ^ )/../-«(-)'^-) (81) 



VI. RESULTS AND DISCUSSION 

First of all we can check that the equilibrium results are the expected ones. It is easy to prove that, independent 
of the dynamics, the effective temperatures are just the temperature of the bath: 

Tls = 2Ei -eQ = T, T^^ = ^^2 -eQ = T. (82) 

Because in equilibrium the energies of the subsystems are the same (see (pO|)). 
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A. Sequential case 



In the ofF-equilibrium case the results are more interesting. It is easy to verify the foUowing expressions for ( |l8| ) up 
to order e^ : 



Tl^ = 2E^{s) 



y4^^ - eQ{s) + (eQis) - .^f^)^ T dt7H.(0e/''^^-^^^-^-^^»'^^ (83) 

JrAs) OS i^iK2' IrAs) Js 



T^^^2E2{s) + 



dE2{s) 



eQ{s) 



fR.{s) ds ~-^'^' ' '~^-^'^' - K,K2 

The expression for the total effective temperature is just 



(eQis) - e^ ^) ^ /* dt'fn, it')ei:' (/«^ (^)-/«. (^))'^^ 



/i?2 (S) 



-'off — 



as 



dCiit,s) 
ds 



Giit,s) + G2{t,s) 



(84) 



(85) 



where the correlations are given by expressions ( pl| , [4^ ) and the responses are given by ( p^ , p9| ). 

From the equations (83|8J) we can see immediately that the effective temperatures are well defined in the regime 
in which the ratio | is finite. Otherwise, the last term in the right hand side of ( |83| ) and (Q) would diverge. At 
zero temperature and up to order e^ in the coupling constant, we have found that Ei, E2,Q decrease logarithmically 
implying that both fn-^ and Jr^ decay like 1/t. Now let us consider t, s both large but t — s ^ s. For a weak coupling 
(i.e e « 0) the value of the effective temperatures are, in the limit s — > cx) but with - finite: 



Tiff 



2Eiis)+0{e^) 



KfK^Al 



8{KiK2-e^)log{t) 



(86) 



T, 



off 



2E2{s) + Oie^) 



KiKlAl 



8{KiK2 - e^)log{t) 



(87) 



This yields in the s —> oo limit a non vanishing relative difference T^g/T^g — 1. This is a consequence of the 
fact that the two energies are different in the long-time regime. Note that each effective temperature verifies the 
equipartition theorem in the limit of long times as expected. The physical interpretation is clear: each system 
is relaxing towards its equilibrium state slowly and at any time we can consider that the systems are at "quasi- 
equilibrium" at their corresponding effective temperatures. Obviously the concept of "quasi-equilibrium" is meaningful 
in a time window smaller than the characteristic time-scale in which the system relaxes (i.e. during this time-scale 
the effective temperatures do not change), hence we need to impose that t/s is finite. 

Let's think now about the global system. As we have seen, the energies for the two systems remain different even at 
infinite times. This can be explicitly seen in figure (m where we show how the relative difference between the energies 
(or the effective temperatures according to (86]87|)) increases monotonically as a function of time (for late times) for 
any value of e. We may then conclude that a coupling in the Hamiltonian is not enough to reach an equalization of 
effective temperatures. 




FIG. 4. Relative energy difference E1/E2 — 1 for two systems with Ki = 2,7^2 ~ 1, Ai — 1, A2 — 1 and different values of 
e = 0.1, 0.3, 0.5, 0.7, 0.9. Note that the relative difference increases with time. 
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This difference of tlie two effective temperatures implies tliat tfiere are some degrees of freedom hotter than others. 
One can then imagine that there is always some kind of heat transfer or current flow going from the "hot degrees" of 
freedom to the "cold" ones. Then, one may ask why the effective temperatures do not asymptotically equalize. The 
reason is that the off-equilibrium conductivity may vanish with time fast enough for the heat transfer not to be able 
to compensate such difference. In this situation, if we now compute the total effective temperature (p3) for the whole 
system we see that in the off-equilibrium regime this temperature does not coincide with the sum of the energies 
of the systems. This fact fortifies the definition of the effective temperature using the FDR (|l|) in off-equilibrium 
systems. For two systems in "local" equilibrium at two different temperatures, despite the fact that each system 
verifies FDT, the sum of the two systems never verifies FDT unless the two temperatures are equal. In our case, we 
have two systems which are in "quasi-equilibrium" at two different effective temperatures, so the T^g would never be 
the sum of the two energies unless the two effective temperatures T^g, T^g were the same. In other words, two systems 
thermodynamically stable at different temperatures are not globally stable when put in contact. 

B. Parallel dynamics 



The effective temperatures (18) can be exactly computed to order e^. In the equilibrium regime, both the full 
expression derived from (151 15(;.77,78) and the general approximate solutions (pq,pq,pO^,plf) up to order e^ computed 



in the equilibrium regime yield the bath temperature for the three effective temperatures. In the non-equilibrium 
case, up to order e^, the results are: 

The expression for the total effective temperature is just: 

dCi{t,s) dCi{t,s) 

Gi(i, s) + G2[t,s) 

where the correlations are given by expressions ( |65| , |66| ) , and the responses are given by ( pO| , ^ . 

As in the case without coupling, the interesting dynamics is when the temperature of the bath is zero. In this 
case, the energies and the overlap decay to zero logarithmically which implies that /B„(t) vanishes like l/t. A careful 
evaluation of the integrals contributing to the e^ term shown in equations (p8,p3) reveals that they are a function of 
t/s which stays finite provided that ratio is finite. As we discussed in the previous uncoupled or sequential case the 
effective temperatures (|8[|9|) have full sense when we consider times t/s finite so no appreciable transfer of energy 
between the two systems has still occurred. 

It is clear from the asymptotic expressions for the energies and the overlap that in the long-time limit (s — > oo) : 

While the energies themselves decay as l/log(s) the relative difference {Q{s)/Ei{s) — £ j^^j^ ) decays like l/log(s). 
Up to order e^ we may write, in the limit s — > oo (with - finite): 

T^'g = 2Eiis) ~ eQis) , T^g ^ 2E2{s) - eQ{s), (92) 

because the asymptotic values of the £'i(s) and E2{s) are the same the effective temperatures for the subsystems 
become identical in the long-time limit. Note that the case with dynamic coupling or parallel dynamics is qualitatively 
different from the case without dynamic coupling or sequential, because now all the degrees of freedom are at the 
same effective temperature in the long-time limit. Moreover, if we consider the global system it is easy to prove that 
the total effective temperature defined in ( |l8| ) is, in the limit s — > cxd with - finite: 

Ti^ff = T,'g = T^g = 2E{s) - eQ{s) (93) 
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where E = Ei = E2 and Q are given by (^6[). This is a consequence of the fact that the energies of the two systems 
equahze due to the dynamic coupUng. Then, the whole system has the same effective temperature and we can define 
an effective temperature for the global system using FDT. The situation is the same as in equilibrium systems. If we 
have two systems in equilibrium at a certain temperature T, FDT not only holds for each subsystem but also holds 
for the whole system bringing the temperature of the bath T. At higher-orders in e we expect that all terms with be 
subleading for t/s to be finite and asymptotically all three temperatures coincide. 

If we restrict to the case in which the coupling constant vanishes, e = 0, then the systems are still coupled only 
through the dynamics and we obtain the same qualitatively results: 

T<fff = 7;iff = T^f = 2E{s) (94) 

with E{s) « sio (s) • ^^ conclude that the dynamic coupling does not allow the presence of more than one effective 
temperature in the whole system because even in the absence of explicit coupling in the Hamiltonian, the dynamics 
itself makes the energies to equalize in the long-time limit regime. 

VII. CONCLUSIONS 

In this paper we have solved exactly the dynamics of two systems of harmonic oscillators. We focused our attention 
on the concept of the effective temperature defined through the FDR eq.(|l|). The effective temperature, a parameter 
defined by a relation of the correlation and response functions, has been introduced in the context of glass theory in 
order to understand the physics behind the dynamic behavior of these out-off-cquilibrium systems. In this paper we 
hope to have clarified some aspects behind the physical meaning of this effective temperature. 

We have studied two types of couplings between the two subsystems of oscillators, both in an aging state, finding 
that the way we couple them is crucial for the validity of the zero-temperature law in the off-equilibrium regime to 
hold. The two cases we studied are the dynamically uncoupled or sequential case and the dynamically coupled or 
parallel case. In short, for the sequential case the coupling between the variables of the two subsystems in the resulting 
dynamics arises only through the Hamiltonian term eQ. For the parallel case, the variables of the two subsystems are 
simultaneously updated leading to further interaction between the two subsystems (on top of the eQ coupling term 
in the energy). 

We have discovered that for the dynamically uncoupled or sequential case the two subsystems asymptotically 
reach different effective temperatures which never equalize. So the whole system is divided in two parts, each part 
characterized by its own effective temperature. The explanation for this odd behavior lies behind the time dependence 
of the off-equilibrium thermal conductivity which decays very quickly to allow for an asymptotic equalization of the 
two effective temperatures. This raises the question whether different interacting degrees of freedom do eventually 
reach the same effective temperature in the asymptotic regime, condition tightly related to the validity of the zeroth 
law for the off-equilibrium aging state. Our conclusion is that the zeroth law is probably valid but hardly effective 
due to the very small conductivity between the twQ subsystems in the aging state. A calculation of the thermal 
conductivity in this model will be shown elsewheretj and reveals that it decreases very quickly with time, the heat 
transfer being unable to compensate for the difference of the effective temperatures of the two subsystems. 

For the dynamically coupled or parallel case, the two effective temperatures equalize and the two subsystems are 
in a sort of thermal equilibrium between them in the aging state. Consequently, the union of the two subsystems has 
an effective temperature which coincides with the temperature of each subsystem. In this case, the direct coupling 
of the two subsystems through the parallel dynamics makes the conductivity much larger than in the sequential case 
so in this case a zero-th law for the aging state is effective and holds. In fact, these results are also valid when we 
consider the particular case e = in which the dynamics in itself is enough to equalize the effective temperatures. 

From these two type of couplings the first one is the only realistic. Dynamics in real structural glasses involves 
short scale motions of atoms and coupling between the different degrees of freedom occurs at the level of the energy 
or Hamiltonian and never at the level of the dynamics (at least in the classical regime). The results of this paper 
explain then why different degrees of freedom in structural glasses can stay at different effective temperatures forever. 
The off-equilibrium conductivity or heat transfer between the different degrees of freedom is small enough for the 
equalization of the effective temperatures associated to the different degrees to never occur. This explains why when 
we touch a piece of glass we feel it at the room temperature. The heat transfer coming from the hotter non-thermalized 
degrees of freedom is extremely small. Before finishing we must note one particular feature of our model. All the 
calculations were done at zero temperature where the energy vanishes asymptotically. The fact that the energy (and 
consequently the conductivity) of the system is exhausted in the asymptotic limit can lead to a pathological behavior 
not present in structural glasses at finite temperature. Nevertheless, the fact that the thermal conductivity vanishes 
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much faster than the energy itself, suggests that the vanishing of the conductivity is not related to zero-temperature 
dynamics. 

In the present calculation we have focused on the interaction between two subsystems, both in the aging state. 
When one of the subsystems is in an aging state and the other is in equilibrium the analysis proceeds similarly, the 
conclusion being that the non-thermalized subsystem determines the rate of heat transfer and hence the measurement 
of the effective temperature. The value of the effective temperature measured by a themiometer and other related 
questions can be analyzed in detail in the present model and will be presented elsewheretj. 

To conclude, although a zero-th law for non-equilibrium glassy systems may hold, it is hardly effective because of 
the small energy transfer occurring between degrees of freedom at different effective temperatures. It would be very 
interesting to pursue this investigation further by studying other solvable examples and showing that what we have 
exemplified here is a generally valid for structural glasses as well as for other glassy systems. 
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APPENDIX A: A SHORT REVIEW OF THE BPR MODEL 

In this appendix we show how derive the results for the correlation and the response functions in El in order to 
understand the techniques we will use throughout this paper. In that model, the system is constituted by N uncoupled 
harmonic oscillators which evolve with Monte Carlo dynamics. The energy of this system is 

E{x,) = IkJ2^1 ■ (95) 



The result for the dynamical evolution for the energy is: 

dE aa fl-AEP 



dt 2 \ acP 
where we have defined the quantities 



f{t) + erfc{a)] , (96) 



1 KA'^ 2 f°° 

Qc^-KA'^ , "=Y76^' ^^f^i^)"^—^ exp{-x'^)dx, (97) 

/ KA^ \ 

f{t)=a,(3expi~/3^-{l-2E{t)/3)\erfc{a{t){4E{t)P~l)) . (98) 

The stationary solution is just E = -^T as expected. Another important quantity is the acceptance rate which is 
the number of accepted Monte Carlo movements at a time t: 

^(.)4(M+e./c(«)) . (99) 



In the same way we can compute the equation for the correlation function defined as: 

N 

N 



C(i,s) = l^x,(i)x,(s) , (100) 

1=1 

and the evolution of the correlation function is given by the equation 

dC{t,s) 



dt 



= -f{t)Cit,s) , (101) 



where the quantity f{t) has been previously defined in (pq). The solution for the correlation function (which depends 
explicitly on two times) is: 

C{t,s)^^^exp(^^Jj{x)dx^ , (102) 

where we have to add the initial condition C{s, s) = — f^ ■ In order to compute the equation for the response function 
defined by, 

Oii..)=('-^) , (103) 



we have to consider the Hamiltonian perturbed by a small external field 

i—l i 

Then we compute the dynamical evolution for the magnetization, which in our model is defined by: M = jj X]i=i ^i 
yielding: 
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^^^(') -(M{t)-^]fEit)eit-.s) , (105) 



dt \ " K 

where we have defined a 'new' energy as: 

E' = E-Mh+—- . (106) 

2iv 

The quantity fE'{t) is identically defined as in (^8|) but with the new energy E' . Note that in the case in which 
h = the magnetization will always be zero because of the initial condition we consider, i.e M(0 ) = 0. Also, when we 



compute the evolution for the response function, the first term in the right hand side of the ( 105 ) is just the response. 
Th en w c have to analyze carefully the second term (which is proportional to the external field) in the right hand side 
of ( |l05| ). First of all, we write 

Git s) - ^^^ ^^M ^ Ziz (107) 

^^*''>~ hAs ' dt ~ hAs ■ ^^°^> 

We consider the variation of the magnetization as follows 

AM ^ d fdM(t)\^^ , ^ dM(t) ^^ , , „„, 

^ = A.^(^je(i-.) + A.^^(t-.) , (108) 

and by keeping only the linear term in h in the second term of the r.h.s in (|10§|) we get 



^ = As^'-§^eit-s) + As^fmt-s) . (109) 



The first term in the r.h.s of (109) is obviously zero and only the the second term gives a non-vanishing contribution. 



So the evolution for the response function is 



^^^''"^ =-f{t)(G{t,s)-is{t-s)\ , (110) 



dt ^ V ^ y V w ^ 

whose solution is 

G{t, s) = M cxp (^- ^ f{x)dx^ e(t - .s) . (Ill) 

Now, we are in position to compute the effective temperature based on the violation of FDT 

Note that the effective temperature only depends on the smallest time s. This feature is due to the simplicity of the 
model. In this model, due to the finite amplitude of the Monte Carlo movements the system never reaches the ground 
state {xi — 0}. In fact, Monte Carlo dynamics induces entropic barriers which manifest as activated behavior for the 
relaxation time. The interesting dynamics is found when we study the relaxation of the system at zero temperature. 
To obtain the dynamical equations at zero temperature we have to consider only the negative changes in the energy. 
It can be seen that in the long time limit the relaxation of the energy is logarithmic 

. ^ KA^ 1 , , 

16 log(^)+ilog(log(^)) 

moreover, we obtain the following asymptotic behavior of the function f{t) and the acceptance rate: 



,/ ^iog(iog(^))y 



log(7l) / 4tlog(^) 

For the long time behavior of the correlation and the response functions we obtain to leading order in log(s)/ log(t) 



2E{s)^ 



^^"^'W .. ^^ 1 hi^y 



C[t,S) — — — C„orm(^:S) J L/„orm(^j S) ~ —/ ^ j '''' '^^ ~ 177 7 \"l ' (H^) 

^ tW(^) ^^Vlog(^) 
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APPENDIX B: SOLUTION OF THE DYNAMICALLY UNCOUPLED OR SEQUENTIAL CASE 

In this appendix we show explicitly the detailed calculations for the case in which we sequentially update the two 
subsystems. Note that each subsystem is updated in parallel but no simultaneous updating of the whole system is 
performed so there is no direct coupling of the two subsystems through the dynamics but only through an explicit 
coupling term eQ in the Hamiltonian. We have to take into account this fact when we compute the distribution 
probability for a change in the energy. The Hamiltonian we have to deal with is 
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where Ei and E2 are the energy of the bare systems while Q is the overlap between them. 

The Monte Carlo updating procedure is the following. First all the Xi are moved to Xi + ri/\/N where the r^ 
are random variables Gaussian distributed with zero average and variance A^. The move is accepted according to a 
rule defined by an acceptance probability W{AE) which satisfies detailed balance: W{AE) ~ W{—AE) exp(— /3A_E), 
where AE is the change in the Hamiltonian. Later all the yi are moved to yi + Si/\/N, where the Si are random 
variables Gaussian distributed with zero average and variance A\ . The same transition probability is now applied for 
the yi variables. This sequential updating of the xi and yi variables is then iterated. Note that the coupling in the 
dynamics only appears through the change eAQ of the total energy. 

Now we compute the distribution probability of a change in the energy of the first system. This probability 
distribution can be expressed 
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and in the same way we can compute the probability for the other system: 
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Using the integral representation of the delta function: 
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with the quantities 
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Note that due to the explicit coupling e the probability of a change in the energy of one system not only depends 
on their energy, but also on the energy of the other system and the overlap. Now, we can compute the evolution of 
the energies 
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We compute the equation for the evolution of the overlap in two steps. The first is the change in the overlap when 
the variables of the first system are moved; and the second one is when the variables of the second system are moved. 
So we must to compute two joint probability distributions 
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Then we compute the evolution equation for the overlap in each step and sum the two equations 
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The solution of these equations is 
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which yields the final equation 
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with the quantities defined in (E7J97). In the same way we can compute the equation for the correlation and overlap 
functions defined in ( |ll| , |l2|) . To compute their evolution equations we must evaluate the joint probability of a change 
in the energy and a change in the correlation function. Note that when we consider the change in the variables Xi we 
have to consider the energy of the system one, and when we consider the change in the variables yi we have to take 
into account the energy of the other system. The joint probability can be decomposed into the probability distribution 
for a change in the energy multiplied by a conditional probability 
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Then the evolution for the correlation functions can be computed using 
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For the response functions we have to compute the dynamic evolution equations for the magnetizations. We consider 
an external field coupled to each system, so the new Hamiltonian is 
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We define the magnetizations as follows 
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Then, we perform the same steps as we did for the other quantities. First of all we have to compute the joint 
probability of a change in the magnetization and a change in the energy. For example, for computing the response 
function for the first system we make /12 = and hi ^ 0; then we compute the joint probability distribution for a 
change in Mi and Ei . After that we can obtain the evolution for the magnetization of this system 
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Note that in this case we are considering ft,2 = but still the equation for Mi depends on M2. For the sequential 
updating procedure we have to consider the evolution for M2{t) with /ii = and /12 7^ which is, by symmetry 
considerations 
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We finally get the equations for the four different response functions using the same procedure we followed for the 
single system (see Appendix A). This yields 
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In order to compute the effective temperatures we shall use a perturbative expansion in terms of the coupling 
constant described in Appendix C. 

A. Equilibrium results 

In equilibrium the matrices for correlations and responses can be exactly diagonalised. The results are 
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with the values of the constants 
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and the two eigenvalues: 
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The results for the other two correlation functions have the same form 
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with the new values of the constants 
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and the same expressions as in eqs.(155) for Ai, A 
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APPENDIX C: SOLUTION FOR THE OFF-EQUILIBRIUM CORRELATIONS AND RESPONSES 

IN THE INTERACTION REPRESENTATION 

In general we have to solve the following equation 

Yt = ^W"^ ' (160) 

with the initial condition v{t) — v{s). A(t) is the matrix with the time-dependent coefficients of our problem. It can 
be decomposed as: 

A{t) ^ Aoit) + eAj{t) , (161) 

where ^o (0 i^ the diagonal part and Aj (t) is the interaction part of the matrix. We work in the interaction repre- 
sentation. Therefore we start by doing the transformation 

w{t) = (exp (- / Ao{t')dAj v{t) . (162) 

The derivative of this new vector is simply: 
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Now we must solve (164) with the initial condition w(s) ~ v{s). The formal solution for this equation is 

w{t):^w{s)+€ I Bi{t')w{t') , (166) 

or equivalently 

w{t)^Texp{e I Bi{s'))w{s) . (167) 



Where T stands for the time ordered product. This equation can be iterated and solved to any order in e. Up to 
order e^ we find 

Order zero : w{t) — v{s) , (168) 

Order e: w{t) ^ v{s) + e / Bi{t')v{s) , (169) 

Order e'^ : w{t) ^ v{s) + e f Bi{t')v{s) + e^ [ dt' f dt" Bi{t')Bi(t")v{s) . (170) 



This is the procedure we have used in order to obtain the equations for the responses and correlations for the 
dynamically coupled and uncoupled cases. 
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